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Abstract: The statistical approach for description of molecular cloud substructure, 
proposed by Donkov, Veltchev and Klessen (2011, 2012), allows for alternative 
models, operating with different type of objects: an ensemble of clumps or a larger 
cloudlet. We demonstrate briefly the predictive power of both models, applied to 
molecular emission and dust extinction studies of Galactic clouds. 



1 Physical description of clumps 



Clumpy, dense structures in molecular clouds (MCs) and/or in the cold neutral medium are 
direct progenitors of protostars or of protostellar clusters. While numerous algorithms for 
identification of coherent regions on column- density maps or in the position-velocity space 
have been suggested (e.g. , IStutzki fc Giistenl . Il990l ; IWilliams. de Geus fc Blitzl . Il994l ; 
Mensh'chikov et al.l . |2010| ). the issue whether and how those objects represent the existing 
clumps is open. Therefore it is useful to i nvent a statistical approach for their description 
flDonkov. Veltchev fc Klessenl . l201lL \20lj ) . 



A physically sound basis for such approach is to consider clumps as condensations of 
turbulent (shock) origin formed and shaped during the early MC evolution while gravity 
slowly takes over at small spatial scales. Assuming that compact (quasi-spherical or cubic) 
clumps with sizes /, significantly lower than the scale of their generation L, and of various 
densities p in respect to the local mean density (p) l populate the whole considered volume 
we adopt a natural physical relation for their masses: 



m 
m 



P 



(Ph 



where mo and Iq are units of normalization. The building blocks of the model stem from the 
physical consideration of: i) isothermal, isotropic, fully saturated turbulence in a steady state; 
and, ii) the role of gravity for the energy balance at different spatial scales. 
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2 Basic assumptions of our statistical approach 



1) Turbulent cascade in isothermal medium in the scale range .5 5, L S, 20 p c, 
conditioning scaling laws of velocity and density according to iLarsonl (Il98ll ): 
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2) Lognormal probability density distribution (PDF) at each scale L with stddev 

a 2 = ln(l + b 2 M 2 ) 



where b is the turbule nt forcing parameter (IFederrath. Klessen &: Schmidtl . 12008 



Federrath et al.l . 120101 ) and Ai is the Mach number. 



3) Mass-density power-law relationship: 
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x < 0, x = x(L) 



4) Self-similar density scaling as expected for turbulent structures: 



5) 'Virial-Hke' equipartition of energies: 

\W\~f sk E kin , l</ gk <4, (wkin{l-4}) 
\W\ ~ 2£; kin + 2£th, (wkin2th2) 
|W| ~ 2^ kin + E mag , (wkin2mag) 



3 Alternative models within the chosen approach 

They are built according to the type of object for which the assumptions of mass-density 
relationship (3) and of a chosen equipartition of energies (5) hold: 



Stat istical ensemble of clumps ('ensemble model'; iDonkov. Veltchev fc Klessenl . 12011 



20121 ): Populations (ensembles) of clumps generated at abstract spatial scale L obey a 



mass-density relationship with power-law index x(L), derived from equipartition 
relation for their representative members ('typical clumps'). 



• Cloudlet, defined by a density cut-off ('cloudlet model'): A chosen energy equipartition 
(assumption 5) determines a density threshold in respect to the mean local density (p)l- 
Then the relationship (3) with a power-law index x(L) describes the intrinsic 
mass-density scaling of structures built into the delineated cloudlet. 
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Figure 1: Predictions of the 'ensemble model' applied to Orion B: (left) Model (line) and 
estimates from dust extinction of the comple x subs tructure in terms of x(L); (right) Modelled 
and observational (CO maps; Kramer et all Il998l ) clump mass function (C1MF). The mass 
range corresponds to spatial scales of clump generation 0.7 < L < 10 pc (shaded areas). 



3.1 Clump mass function from the 'ensemble model' 



The assumption of 'virial-like' energy equipartition wkin2 or wkin2th2 yields a 
description of MC structu r e in a greement with dust-absorption maps of large clouds 
flLombardi. Alves fc LadaL l20ld ). 



The modelled high-mass CIMFs exhibit power-law slopes T ~ — 1, as expected for 
fractal clouds (jElmegreenl . 119971 ) . When only gravitationally unstable clumps are 
considered, the slope of their time-weighted high-mass mass function is similar to that 
of the stellar IMF: T 1.3. 



3.2 Clump mass function from the 'cloudlet model' 

• The intrinsic mass-density scaling within the cloudlets is consistent with that of MC 
fragments delineated by the dendrogram technique in some Galactic clouds 



flKauffmann et al!l2010h . 



• Further insights into physics of the dendrogram objects would allow for derivation of the 
C1MF within the 'cloudlet model'. 
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Figure 2: Modelled C1MF from various equipa rtition relations, compa red with the observational 
one, derived for a sample of 9 Galactic MCs ( iTachihara et al.l . |2002| ). 



-0.5 



-1.5 



T=10 K, p=0.33, b=0.55 
T=20 K, 0=0.33, b=0.33 
T=20 K, 0=0.43, b=0.55 




1 

Effective size L c [ pc ] 



10 



10" 



I -0.83 10 

r=-i.oo 




10 J 10* 
M t [M Sun ] 



0.1 



0.01 



Figure 3: Predictions of the 'cloudlet model': (left) Calculated x(L c ) in cloudlets from the 
equipartition wkin2th2, compared with observational r anges for MC fragments in Ophiuchus 
(black) and Perseus (violet), derived from the results of iKauffmann et al.l ( 12010I . see their Fig. 
5); (right) Cloudlet mass functions, derived for the same sets of model parameters. The slopes 
vary from a typical value for CO clumps to that for fractal clouds. 
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